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Abstract 



CO 

The stability of isotropic cosmological solutions in the Bianchi I model is considered. 
We prove that the stability of isotropic solutions in the Bianchi I metric for a positive 
Hubble parameter follows from their stability in the Friedmann-Robertson-Walker 
metric. This result is applied to models inspired by string field theory, which violate 
the null energy condition. Examples of stable isotropic solutions are presented. We 
■ also consider the ^-essence model and analyse the stability of solutions of the form 

> ' 

1 Introduction 

a. 

Field theories which violate the null energy condition (NEC) are of interest for the solution of 
the cosmological singularity problem [TJ [2j [3] and for models of dark energy with the equation 
of state parameter w < — 1 (see [I] [II] and references therein). Generally speaking, models 
that violate the NEC have ghosts, and therefore are unstable and physically unacceptable. 

However, the possibility of the existence of dark energy with w < — 1 on the one hanoS 
and the cosmological singularity problem on the other hand encourage the investigation of 
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models which violate the NEC. It is almost clear that such a possibility can be realized within 
an effective theory, while the fundamental theory should be stable and admit quantization. 
From this point of view the NEC violation might be a property of a model that approximates 
the fundamental theory and describes some particular features of the fundamental theory. 
With the lack of quantum gravity, we can just trust string theory or deal with an effective 
theory admitting the UV completion. 

There have been several attempts to realize these scenarios [18j[T9], 20J. The ghost conden- 
sation model [TH [HJ (221 123] proposed to describe a wide class of cosmological perturbations 
has a ghost in the perturbative vacuum and has no ghost in the ghost condensation phase 
within an effective theory. The new ekpyrotic scenario [201 [2U EHl [26] is a development of 
the ekpyrotic [27] and the cyclic scenarios [21 EE], and it attempts to solve the singularity 
problem, among others, by involving violation of the NEC. Nonlocal cosmological models 



HS1I2SIEDIEI1E21E31E11ESIES1EZIEB1 inspired by the string field theory (SFT) [391 M, E] 



admit a regime with w < —1. 

All of these models possess higher derivatives terms, which produce well-known problems 
with quantum instability [121 S3]- Several attempts to solve these problems have been 
recently performed [HI [26] . A physical idea that could solve the problems is that the 
instabilities do not have enough time to fully develop. A mathematical one is that dangerous 
terms can be treated as corrections valued only at small energies below the physical cut-off. 
This approach implies the possibility to construct a UV completion of the theory, and this 
assumption requires detailed analysis. 

The NEC plays an important role in classical general relativity, in particular, in the 
consideration of black holes and cosmological singularities [HE]. 

The NEC violating models can admit classically stable solutions in the Friedmann- 
Robert son- Walker (FRW) cosmology. In particular, there are classically stable solutions 
for self-interacting ghost models with minimal coupling to gravity. Moreover, there exists 
an attractor behavior (for details about attractor solutions for inhomogeneous cosmological 
models, see |45j) in a class of the phantom cosmological models [IHl H7J 0E]. One can study 
the stability of the FRW metric, specifying a form of fluctuations. It is interesting to know 
whether these solutions are stable under the deformation of the FRW metric to an anisotropic 
one, for example, to the Bianchi I metric. In comparison with general fluctuations we can 
get an explicit form of solutions in the Bianchi I metric, which can probably clarify some 
nontrivial issues of theories with NEC violation. 

Stability of isotropic solutions in the Bianchi models [HI [501 EI] ( see also [22]) has been 
considered in inflationary models (see [531 151] and references therein for details of anisotropic 
slow-roll inflation). Assuming that the energy conditions are satisfied, it has been proved that 
all initially expanding Bianchi models except type IX approach the de Sitter space-time [53] 
(see also [561 EH ESI [59]). The Wald theorem [55] shows that for space-time of Bianchi 
types I- VIII with a positive cosmological constant and matter satisfying the dominant and 
strong energy conditions, solutions which exist globally in the future have certain asymptotic 
properties at t — > oo. It is interesting to consider a similar question in the case of phantom 
cosmology [UJ EDI HE] and string inspired models [THJ [201 EH ESI ED E2] , as well as in the 
case of the ghost condensation models [18J or their modifications [26] . 




The Bianchi universe models [191 [501 151] are spatially homogeneous anisotropic cosmo- 
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logical models. There are strong limits on anisotropic models from observations [63, 64J. 
Anisotropic spatially homogeneous fluctuations have to be strongly suppressed, and models 
developing large anisotropy should be discarded as early or late cosmological models. 

In this paper we consider the stability of isotropic solutions in the Bianchi I metric in the 
presence of phantom scalar fields. There are two classes of models whose stability we analyse 
in this paper. The first class includes the one phantom scalar field models of dark energy, 
which admit exact kink- type or lump- type solutions [601130] . For this class of models we also 
analyse the stability with respect to small fluctuations of the initial value of the cold dark 
matter energy density (compare with [46J). The second class includes models with a scalar 
field 4>, which have exact solutions <fi ~ t, for example, the fc-essence models [651 EEl E7J 168] . 
in particular, ghost condensate models [T8 l I2T1 126] . 

For both classes of models we prove that the stability of the solutions in the Bianchi I 
metric is equivalent to the stability of the corresponding solutions in the FRW metric. The 
stability of a kink or lump solution in the FRW metric means the stability of the fixed point 
that the solution tends to. Using the Lyapunov theorem [691 El] we find conditions under 
which the fixed point and the corresponding kink (or lump) solution are stable. In these 
cases the necessary condition for the exact solution's stability is boundedness of the first 
corrections for the positive time semiaxis. When we can not use this theorem we check the 
boundedness of the first corrections to the exact solutions explicitly. 

The paper is organized as follows. In Section 2 we deal with an arbitrary iV-component 
scalar potential model and a fc-essence model in the Bianchi I metric. We also review the 
Lyapunov theorem and other important statement about stability. In Section 3 we consider 
the stability of solutions which tend to an isolated fixed point in one-field models with the 
cold dark matter (CDM). We find sufficient conditions for the stability of such solutions in 
the FRW and Bianchi I metrics. In Section 4 we find the connection between the first order 
corrections in the FRW and Bianchi I metrics for iV-field models. The corresponding result 
for the fc-essence model is presented in Section 5. In Section 6 we present examples of stable 
isotropic kink and lump solutions in SFT inspired cosmological models. We also analyse the 
first order corrections for solutions, which are proportional to time. In Section 7 we make a 
conclusion and propose directions for further investigations. 

2 Setup 

2.1 The Bianchi I cosmological model with scalar and phantom 
scalar fields and the CDM 

Let us start with a cosmological model with N scalar fields 0i, 02, • • • , <Pn i n the Bianchi I 
metric 



ds 2 



dt 2 + a\(t)dx\ + a 2 (t)dx2 



+ 




(1) 



The action is 



S = d A x\/—g 




lQnG N 



R 



Y 9 ^ d ^ kdu( ^ k 



k=l 



(2) 
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where the potential V is a twice continuously differentiable function, Gn is the Newtonian 
gravitational constant, A is a cosmo logical constant, and C k are nonzero real numbers. The 
sign of Ck defines whether field (f) k is the phantom field (C k < 0) or the ordinary scalar field 
(C fc >0). 

The Einstein equations have the following form: 

H X H 2 + H X H 3 + H 2 H 3 = 8nG N g, (3) 

H 2 + Hi + H 3 + Hi + H 2 H 3 = - 8nG N p, (4) 
Hi + Hi + H 2 + Hi + H t H 2 = - 8nG N p, (5) 
Hi + Hi + H 3 + Hj + HtHs = - 8nG N p, (6) 



where 



N 



fe=i 

A' 



k=l 

H 1 = ^, H 2 = ^ H 3 = ^ (9) 
ai a 2 a 3 

and a dot denotes a time derivative. 

Note that we couple, in a minimal way, pressureless matter (the CDM) with the energy 
density p m to our model. The equation for the CDM energy density is as follows: 

p m = ~{Hi + H 2 + H 3 )p m . (10) 

Introducing ip k = 4> k we obtain from action (j2j) the following equations: 

4>k = i>k, 



j> k = -(Hi + H 2 + H 3 )^ k -—V; k , 
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where V^ k = k = 1, 2, . . . , N. Thus we get the system of 2N + 4 first order differential 
equations and one constraint (j3J). 

It is convenient to express the initial variables dj in terms of new variables a and /3, (we 
use notations from [71]). subject to the following constraint: 

A + & + ft = 0. (12) 

One has the following relations 

fli (t) = a(t)e ftW , hence, a(t) = (ai(t)a 2 (t)a 3 (t)) 1/3 , (13) 

Hi = H + Pi, and H = ^(Hi + H 2 + H 3 ), (14) 
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where H = a/a. To obtain (1141) we have used the following consequence of (1T2"|) : 

k + k + $3 = 0- (15) 

Note that /3j are not components of a vector and, therefore, are not subjected to the 
Einstein summation rule. In the case of the FRW metric all /3{ are equal to zero and H is 
the Hubble parameter. Following [501 [71] (see also [52]) we introduce the shear 

o* = % + % + (16) 

It is useful to write equations (EI) (J6]) , i fTOl) and ( fill in terms of new variables. 
Using relation (1151) we can write equation ([3]) as follows 

3ff 2 - X -o 2 = 8nG N g. (17) 

Summing equations (SD (JH]) one can obtain 

2if + 3ff 2 + \o 2 = - 8nG N p. (18) 

Therefore 

if + 3H 2 = 4irG N (g-p). (19) 
Note that equations ( fTUl) and ( fill in new variables, 



fc = V*, ^ = - 3#V* - 7^ fe , (20) 



p m = - 3Hp m , (21) 

as well as equation (fl~9]) . look like the corresponding equations in the FRW metric. 
Subtracting (jlj) from (j5J we obtain 

H x + ff x 2 - ff 3 - fff + ff 2 (ff x - ff 3 ) = 0. (22) 

In terms of if and $ equation (1221) takes the form 

A + 3if /3 a = /3 3 + 3ff /3 3 . (23) 

Using (1231) and (jT5l) we obtain the following equations 

A= -3i?A, (24) 

| (a 2 ) = - ma 2 . (25) 

Functions if (£) and cr 2 (t) together with 0fc(t) and p m (t) can be obtained from equations 
(TT51)-(|2"T1) and (1251) . If if (t) is known, then /3 4 can be trivially obtained from (!24j) . We show 
in the next section that functions H(t), $i(t), and cr 2 (t) are very suitable to analyse the 
stability of isotropic solutions in the Bianchi I metric. 
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2.2 /c-essence model in the Bianchi I metric 

Let us consider the fc-essence cosmological model, which is described by the action 

S = / A ^(l6^- P ^- A )' < 26) 

where 

X= -gTd^d^. (27) 
The pressure V($, X) is of the form [TH [20] 

P(*,X) = \{p q m - QM) + \W*) + Q 9 {*))X + \m\$)(X - l) 2 . (28) 
Here p q {<&), Qq{&), and M 4 ($) are arbitrary functions of $. The energy density is 

X) = (p,($) + Q q ($))X + 2M 4 ($)(X 2 - X) - P($, X). (29) 

In the Bianchi I metric for $, depending only on time, we have X = $ 2 . The Einstein 
equations are 

H 1 H 2 + H 1 H 3 + H 2 H 3 = 8nG N (£ + A), (30) 
H 2 + H% + H 3 + H% + H 2 H 3 = -8irG N {V-A), (31) 
H 1 +Hf + H 2 + Hl + H 1 H 2 = - 8nG N (V - A), (32) 
H^Hl + H^ + Hl + H^ = -8irG N (V-A). (33) 

From action (1261) we also obtain the second order differential equation for the fc-essence 
field $, which represents a consequence of system (I3"U|) - (I3"3"|) . Indeed, we differentiate ( 1301 
with respect to t and obtain 

(H 2 + # 3 )#i + (Hi + H 3 )H 2 + (#! + ij 2 )# 3 = SttG^. (34) 
Using (1301) - (1331) to exclude Hi, we transform this equation into the following form: 

8 = -(H 1 + H 2 + H 3 )(S + P). (35) 

Substituting explicit forms of £ and V , we obtain 



(2T q + M 4 (3$ 2 - 1)) $ = - T^$ 2 - U; - 2M 3 M' (3$ 4 - 2$ 2 - l) - 

- 2(#i + # 2 + # 3 )$ (T q + 2M 4 ($ 2 - 1)* 



(36) 



where a prime denotes a derivative with respect to $, 

K($) = -*,(*)), T„($) = i(p 9 (*) + (37) 
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The fc-essence model has one important property. For any real differentiable function 
H (t), there exist such real differentiable functions g q (3>) andp g ($) that the functions H,i{t) = 
H (t) and —t solve system fl30|) - fl33|) and, therefore, equation fl36|) . Indeed at = t, 

£ = Q q ($) = Q q (t), V = p,($) = p q (t). (38) 

So, one can obtain from fl30l) - (l33l) 

qM = s^-^oC*) - A ' p«(*) = - eS) - Z^-#(*)- (39) 

Substituting the obtained g q (<&) and in (T28~j) . we see that the system (!30l) - (!33l has a 

particular solution ifj(t) = Ho(t) and = £. 

Bianchi-type models I- VIII coupled to fc-essence matter representing dark energy and 
other matter which satisfies the strong and dominant energy conditions have been considered 
in [72]. A general criterion for isotropization of these models has been derived |72j. In this 
paper we do not assume that the energy conditions are satisfied when considering fc-essence 
models in the Bianchi I metric. 



2.3 A few known facts about stability 

Let us remember a few facts about the stability [69l [70j [73] of solutions for a general system 
of the first order autonomic equations 

y k = F k (y), k = l,2,...,N. (40) 

By definition a solution (a trajectory) yo(t) is attractive (stable) if 

\\y{t) -y (t)\\ ^0 at t -> oo (41) 

for all solutions y(t) that start close enough to yo(t). 

If all solutions of the dynamical system that start out near a fixed (equilibrium) point 

Vh 

F k (y f )=0, fc = l,2,...,JV (42) 

stay near yf forever, then yf is a Lyapunov stable point. If all solutions that start out near 
the equilibrium point yf converge to yf, then the fixed point yj is an asymptotically stable 
one. Asymptotic stability of fixed point means that solutions that start close enough to the 
equilibrium not only remain close enough but also eventually converge to the equilibrium. A 
solution yo(t) of (j4"0]) . which tends to the fixed point yf, is attractive if and only if the point 
yf is asymptotically stable. 

The Lyapunov theorem [69, 70J states that to prove the stability of fixed point y/ of non- 
linear system (|4"0"|) it is sufficient to prove the stability of this fixed point for the corresponding 
linearized system 

y = Ay, A ik = — \ y= y f . (43) 
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The stability of the linear system means that real parts of all solutions of the characteristic 
equation 

'OF 

det ' ~dy~ ~ XI ) \ y=yf = ° 



are negative. 

In the case of a hyperbolic fixed point, i.e. the case when the Jacobian matrix of F at 
the fixed point does not have eigenvalues with zero real parts, one can use the Hartman- 
Grobman theorem [T4"l [75| [76] . This theorem reduces the study of the system of the first 
order nonlinear equations near the hyperbolic fixed point to the study of the behavior of its 
linearization near the origin. 

The case with pure imaginary eigenvalues of the Jacobian matrix of F at the fixed point 
requires a more specific treatment [76J. 

3 Stability of isolated fixed points and kink-type solu- 
tions in one-field models with the CDM 

Let us consider the gravitational model with one scalar field <fi and an arbitrary potential 
V((p), described by action (J2J) at AT = 1. Equations ffT9l) and (1201) for one- field models are as 
follows 

H= - 3if 2 + 8ttG n (V((J)) + A), 

4> = "0, (44) 

i> = - 3fftf - 

This system of three first order equations is valid in the Bianchi I metric as well as in the 
FRW one. Different initial values of a 2 in ffTTj) specify these different cases. 
Let us define 

I = —^—H 2 - ^ 2 - V(<f>) - A. (45) 

From system (|44l) it is follows that the function / should be a solution of the following 
equation: 

i= -6HL (46) 

If the case H(t) = is excluded, then / is an integral of motion of fT4*4|) if and only if 
1 = 0. From (fTTj) we see that 

1 = ttV^ 2 ' ^ 

IoitGn 

so, I is an integral of motion only at a 2 = 0, i.e. in the FRW metric. From ([4*6!) and fT4T|) it 
follows that equation (|25|) is a consequence of (jHJ). 

We are interested in the stability of kink and lump solutions, namely, we consider such 
solutions in which the Hubble parameter tends to a finite value at t — > +oo. In this case 
4>(t) tends to a finite value as well. Thus, there exists a fixed point yf = {H^, ipj), which 
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corresponds to t = +00. We consider the stability of isotropic solutions only, so cr| = and 
$i f =0. It is easy to see that 

V/ = 0, V;(<f> f ) = 0, H 2 f = ~7rG N (A + V(<f> f )). (48) 

To analyse the stability of yf, we present solutions as follows: 

H = Hf + eh(t) + 0(e 2 ) (49a) 

<p = <p f + e(f(t) + 0{e 2 ) (49b) 
ij = e X (t) + 0(e 2 ), (49c) 
$ t =eQ(t) + 0(e 2 ), (49d) 
where e is a small parameter. To first order in e we obtain the following system of equations: 

h(t) = -6H f h(t), (50a) 
<p(t) = X (t), (50b) 

X(t) = - 3H fX (t) - ^Vftfffr. (50c) 
Equation fl50al) has the solution 

h(t) = b e- 6H t\ (51) 

where &o is a constant. 

From (150bl) - (l50cl) we obtain the following solutions: 

9(7 

. at VJ (0,) ^ and (<P f ) ? —H 2 f , 



4 

v? (t) = e- 3 ^ t / 2 ( J D 1 + J D 2 t), (53) 



at y; (<f> f ) = 0, 



= Di - ^fD 2 e- 3H f\ (54) 



where D\, D\, and Z) 2 are arbitrary constants. 

Using the Lyapunov theorem we state that fixed point yf is asymptotically stable and, 
therefore, the exact kink-type or lump-type solution yo(t) is stable if 

V" (6 f ) 

* X > and H f > 0. (55) 
C 

Namely, y/ is 
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• a stable focus at > 

• a stable node at §#| > V *£" > 0, 

• a stable improper node at VI = ^j-Hj. 
From ( J24l) we obtain 

Ci(t) = Qe- 311 ^, (56) 

where Cj are constants. If Hf > 0, then Q and a 2 tend to zero at t — > oo. Thus the 
obtained conditions (1531) are sufficient to prove the stability of isotropic fixed points both in 
the Bianchi I and in FRW metrics. 

At V'l (0/) = or Hf = we need an additional analysis of stability, because the 
Lyapunov theorem does not state the correspondence of the behavior of solutions to the 
initial system (1441) and the obtained linear system (1501) . 

At Hf < the fixed point yf is unstable. Note that both hit) and d(t), as well as <p(t), 
tend to infinity at Hf < 0. 

Let us introduce the CDM into our model with a scalar field. Adding to system (jUj) the 
CDM energy density p m and the corresponding equation (121]) . we get 

H= -3H 2 + AtxGn (2V{(f>) + 2A + p m ) , 

</> = V, 

1 , 57 

Pm = ^>H p m . 



Let us consider the possible fixed points of system (joTj) . From the last equation of this 

system, it follows that at the fixed point we have either Hf = or p m f = 0. Substituting 
(pED and 

Pm(t) = Pmf + ep m (t) + 0(e 2 ), (58) 
into the system (15T|) . we obtain the following system in first order to e: 

Pm{t) = - 3H f p m (t) - 3p rnf h(t), (59a) 

h{t) = - QH f h(t) + 8nG N p m (t), (59b) 

<p(t)= X (t), (59c) 

X(t) = - 3H fX (t) - iy;(0/)<p. (59d) 



It is easy to see that the third and fourth equations of ( 1591) coincide with the corresponding 
equations of system (jSJ). Therefore, the case Hf = can not be analysed by the Lyapunov 
theorem. Let us prove that condition (|55|) is sufficient for the stability of fixed points for 
models with the CDM. First of all, from Hf ^ it follows that p m f = 0. Solving the first 
and second equations of (1591 . we obtain 

p m (t) = &ie~ 3 ^, h(t) = boe^ + ^-e-^\ (60) 

on f 
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where b\ is an arbitrary constant. 

We come to the conclusion that if conditions fl55l) are satisfied, then the solution, which 
is stable in the model without the CDM, is stable with respect to the CDM energy density 
fluctuations as well. 



4 Connections between the first order corrections to 
isotropic solutions in the FRW and Bianchi I metrics 

In the previous section we studied one-field models and the first corrections near a fixed 
point. In this section we consider the first corrections of an arbitrary isotropic solution. 

We consider an iV-field cosmological model, which is described by action (T5]) and the 
Einstein equations P|)- (fTTl) . In this section we do not assume that the isotropic solution 
tends to a fixed point. We do not prove the stability of solutions, we only analyse the 
first corrections in the FRW and Bianchi I metrics. To apply the Lyapunov theorem it 
was convenient to consider functions H and f3{ instead of Hi. In this section we return to 
functions H^. 

To study the stability of this solution, we present solutions whose initial conditions are 
close to the isotropic one, in the following form: 

Hi(t) = H {t)+eh t {t) + O{e 2 ), (61) 

Mt) = 4>v k {t) + ey k {t) + 0{e 2 ), (62) 

M*) = ^0k(t)+e X k(t) + O(E 2 ), (63) 

Pm(t) = Pm0 (t)+ep m (t) + O(e 2 ), (64) 

where % = 1, 2, 3 and k = 1, . . . , N. From (!3j)- (|TT|) we obtain to zero order in e the system of 
Einstein equations and equations of motion in the FRW metric. To first order in e we have 
the following system: 

0fc = Xk, (65) 



1 

Xk = ~ {hi + h 2 + h 3 )?p 0k - 3H Xk ~ tt^2 V 'L<$> m (0o) Vm, (66) 

k m=l 

Pm = - {hi + h 2 + h 3 )p mQ - 3H p m , (67) 



h 1 + h 2 = -3H (hi + h 2 ) + 87rG N J2( V k( ( l ) o)'Pk-C k ^okXk) , (68) 

fe=i 

TV 

hi + fo = -SH {h 1 + h 3 ) + 87rG N J2( V k( ( l ) o)'Pk-C k ^ ok x k ) , (69) 

fe=i 

TV 

h 2 + h 3 = -3H (h 2 + h 3 ) + 87rG N J2(K k ^^-C k ^ 0kXk ) ■ (70) 



k=l 
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From equations (l6"5j) - (IT0j) we get 



k(t) - h 2 (t) + 3# (t)(M*) - h 2 {t)) = 0, (71) 
h^t) - h 3 (t) + 3# (*)(M*) - h(t)) = 0, (72) 



and we also have 

N 

Hoih + h 2 + h 3 ) = AttG n J2 (Ck<PoWk + Vj fc (0 O ) Vfc) • (73) 

k=l 

Theorem 1 

oo 

Let Ho(t) be a smooth function bounded at all finite values of time and J Ho(r)dr be 

o 

bounded from below, in other words, this integral is equal to either a finite number or plus 
infinity. Functions hi(t), h 2 (t), h 3 (t), p m (t), and ipk(t), which are solutions of (65\) -(7(fy, 
are bounded if and only if isotropic solutions, namely, solutions, which satisfy the condition 
h\(t) = h 2 (t) = h 3 (t), are bounded. 

Proof. It is trivial that if the full set of solutions includes only boundary functions, then 
any subset which satisfies an additional condition includes only boundary functions. Let us 
prove that the boundedness of isotropic solutions is not only a necessary condition, but also 
a sufficient one. 

From equations (ITT]) and ( 1721) we obtain: 

t t 

-3 [ HoMdr -3 f HoMdr 

h 1 (t)-h 2 (t) = (h 1 (0)-h 2 (0))e o , h 1 (t)-h 3 {t) = (h 1 (0)-h 3 (0))e ° .(74) 

t 

So we obtain that if the integral J H (r)dT is uniformly bounded from below, then 

o 

anisotropy is bounded at all t. Note that in the most of cosmological models H (t) > for 
all t > and the anisotropy tends to zero at t — > oo. 

Using (1741 . one can express h 2 (t) and h 3 (t) via hi(t) and reduce system (!68|) - (!70|) to one 
equation. System (!65|) - (!70|) takes the following form: 

/ -3j H (r)dr\ f* . \ 

2H I 3h 1 - C e o \ = 8 nG N f£ C k <j> Qh <p k + ^ (0 O ) <Pk J , (75) 

/ N N \ 

2/ii + QHohx = 8-kGn ^ (0 O ) <fk ~ ^ c k4>ok<Pk , (76) 



v fc=i fc=i 



where C = 2/i x (0) - h 2 (0) - h 3 (0). 
Let us introduce a new function, 



C n ~3fH (r)dr 

h (t) = hit) - -fe o . (77) 
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It is easy to check that 

3h (t) = h 1 (t) + h 2 (t) + h 3 (t). (78) 

System (175j)-(|76l) in terms of ho and ipk coincides with the system of equations ( I65i) -( l70l) 
with hi(t) = h 2 (t) = h 3 {t) = h (t). In other words, we obtain that the functions (fk(t) in the 
Bianchi I and FRW metrics are the same. Functions hi(t), h 2 (t), and h 3 (t) differ from the 
correction for the Hubble parameter ho(t) on a finite value. Thus the theorem is proven. 

Note that Theorem 1 connects the stability properties of the FRW and Bianchi I metrics 
not only for solutions which tend to a fixed point, but also for solutions which tend to infinity 
at t — > oo. Examples of such solutions in the cosmological models are presented in Sections 
5 and 6. 



5 Stability of solutions in the /c-essence model in the 
Bianchi I metric 

5.1 First order corrections 

Let us consider the first order corrections in the fc-essence model. Substituting 

£ = £ Q + e£ 1 + 0(e 2 ) } p = V + eV 1 + 0{s 2 ), (79) 

in (l3"0]) - (13"3l and expanding (fBTl) to first order in e, we obtain the following system: 

2H (h 1 + h 2 + hz) = &kG n E x , (80) 

hx + h 2 = -3Ho(hi + h2)-&nG N Vi, (81) 

hx + ha = -3H (h 1 + h 3 )-S-KG N V 1 , (82) 

h 2 + h 3 = -3H {h 2 + h 3 )-87rG N V 1 . (83) 



It is easy to see that equations (ITT!) and (172!) can be obtained from ([8~Ti) - (l53l . therefore, 
formula (1711) is valid for solutions of system (jS0|) - (lS3l . So, it is useful to introduce h by 
formula (I77|) . Because of (1781) we obtain that system (180]) - (1831 in terms of h , E\, and 
V\ coincides with the corresponding equations in the FRW metric, so if Hq satisfies the 
conditions of Theorem 1, then solutions of (180]) - (1831 are bounded if and only if isotropic 
solutions, namely, solutions which satisfy the condition h±(t) = h 2 (t) = h 3 (t), are bounded. 
This means that it is sufficient to calculate the first order corrections for the given background 
solutions in the FRW metric to describe their behavior in the Bianchi I metric. 

5.2 Example 

Let us consider the following example: 

^($) = C 1 + J B$ 2 , p,($) = C 2 - 5$ 2 , M($) = M , (84) 

where B, C\, C 2 , and M are constants. The Friedmann equations are 

3H 2 = 8ttG n (S + A), (85) 
H = -4irG N (S + V). (86) 
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One can check that the following exact solution exists: 



$ (t) = t, H (t) = - AttG n (Ci + C 2 )t, (87) 

if 

C x = - A, B = &kG n (C x + C 2 f. (88) 
Let us analyse the stability of the exact solution, 

$ = %(t)+eV(t), H = H (t) + eh(t). (89) 
The equations for the first order fluctuations, 

h{t) = -87rG N (C l + C 2 + 2M 4 )4>(t), (90) 

9(t) = - 12irG N (^^4^4) K g i + C *M*) + H h (t)} , (91) 
have the following general solution: 

h(t) = d ie 67rGNiCl+C2)t2 + d 2} (92) 

1 I rf l D 67rG N (C 1 +C 2 )t 2 , 2d 2 



87tGat \C 2 + Gi + ZMq L-i + G 2 / 

where d± and c?2 are arbitrary numbers. 

If C\ + C 2 < 0, then H > at t > and the exact solution is stable in the sense that the 
first corrections are bounded functions. Similar solutions, obtained from the SFT inspired 
model, are considered in Subsection 6.5. 



6 Examples of isotropic stable solutions in the SFT 
inspired models 

6.1 String field theory inspired cosmological models 

An interest in cosmological models coming from open string field theories [19] is caused by 
a possibility to get solutions rolling from a perturbative vacuum to the true one. When all 
other massive fields are integrated out by means of equations of motion, the open string 
tachyon acquires a nontrivial potential with a nonperturbative minimum. For the open 
fermionic NSR string with the GSO— sector [ID] in a reasonable approximation, one gets 
the Mexican hat potential for the tachyon field (see [41] for a review). Rolling of the tachyon 
from the unstable perturbative extremum towards this minimum describes, according to the 
Sen conjecture [JT], the transition of an unstable D-brane to a true vacuum. In fact one gets 
a nonlocal potential with a string scale as a parameter of nonlocality. After a suitable field 
redefinition the potential becomes local, meanwhile, the kinetic term becomes nonlocal. This 
nonstandard kinetic term has a so-called phantomlike behavior and can be approximated by 
a phantom kinetic term. Rolling solutions are particular examples of kink-type solutions. 
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It is also interesting to study lump-type solutions, which in particular have no time 
singularity. In was advocated in [301 ED] that such solutions are also available in the SFT 
inspired models. 

In this section we consider the stability of the kink-type and lump-type solutions for the 
SFT inspired cosmological models [SQl EQl EH] under perturbations in the Bianchi I metric. 

In [60l [30] we have considered the SFT inspired phantom models with high degree poly- 
nomial potentials. We consider the stability of the obtained exact solutions in the next two 
subsections. 

In [36] we have considered a nonlocal cosmological model with quadratic potential and 
obtained that exact solutions of this model are solutions of local models with quadratic 
or zero potential. In Subsections 6.4 and 6.5 we analyse the stability of these solutions in 
massless and massive cases correspondingly. The exact solutions in the massive case are 
similar to solutions (187]) in the fc-essence model. 

In the examples we use a dimensionless parameter m 2 ~ M 2 = 1/(8ttGn)- The coefficient 
of proportionality arises when we construct effective cosmological models from the original 
SFT action (see [601 EQ], [36] for details). For convenience, we write the Einstein equations 
for the SFT inspired cosmological models in the following form: 

3 rr2 1 fC^ 2 



We also have 



H = -2 HA -2^ p [-t- V ^- A 
= 4>, (94) 

i> = -3^-^(0)- 

3m 2 p H 2 - - V(cf>) = A. (95) 



6.2 Model with a kink solution and the sixth degree potential 

An exact solution to the Friedmann equations with a string inspired phantom scalar matter 
field has been constructed in [60] (see also [IS])- The notable features of the model are a 
phantom sign of the kinetic term (C = —1) and a special polynomial form of the effective 
tachyon potential: 

W) = ^(l-0 2 ) 2 + T ^ ? 2 (3-0 2 ) 2 . (96) 

Note that this potential has been used in the string gas cosmology [8] . 
System fl9~4|) has the following exact kink-type solution 



4> (t) = tanh(t), H (t) = — ^ tanh(t) (1 - - tanh(t) 2 ) . (97) 

2m 2 V 3 J 

Let us analyse the stability of this solution. At t — > oo solution ( |97j) tends to a fixed 
point, 

H, = -L, *, = 1. (98) 
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It is easy to see that 

Vl(l) = 0, Vft(l) = 2 ( 2 - -L) . (99) 

Using fl55|) . we obtain that solution fl97j) is attractive in the Bianchi I metric at m 2 < 1/2. 
Note that this solution is stable with respect to small fluctuations of the initial value of the 
CDM energy density as well. 

In [nO] we have showed that the first corrections <p(t) and h(t) satisfy the following system: 



h = — jr(l - tanh(t) 2 )(p, 
rrip \ J 



(3 - 4m 2 + 4(m 2 - 1) tanh(t) 2 + tanh(t) 4 ) tanh(t) 

^ = 2m 2 (l -tanh(t) 2 ) ^ ~ ( ^ 

(3 - tanh(t) 2 ) tanh(t) 



1 - tanh(t) 2 
and have the following explicit form: 

<p(t) = 2m\C x (l - tanh(t) 2 ) + 



~ 2 Jit) + (cosh(2t) - l)(cosh(t)) ™i e \^(™^)+v 
+ cosh(2t) + l ' ( 101 ) 



h(t) = G x (l - tanh(t) 2 ) 



2 4m 2 C 2 J(t) 



(cosh(2t) + l) 2 ' 

where Ci and C2 are arbitrary constants, 

rt 2 1 

J(t) = / sinh (r) (coshfY)) 1 - 17 "^ (2 (2m 2 , - l) cosh(r) 2 - l) e 4m2 p c ° shM ' 2 dr. 
Jo 

It is easy to see that if m 2 > 1/2 then at C2 7^ the function ip(t) tends to infinity as 
t — > 00 and, therefore, solution (1971) is not stable. At m 2 = 1/2 we obtain from (110 II) that 

h(t) = (tanh(t) 2 - l) 2 (d - C 2 J 2 ) , 

= - (tanh(t) 2 - 1) (d - C 2 J 2 ) - Ic 2e - tanh(t)2/2 , 

where J2 = J^g-tanM 7 ") 2 / 2 tanh(r)c?r. Thus, <p(t) and are bounded functions at m 2 = 1/2. 
The functions hj have the form 

tanh^ (£) / 2 \ 

^(i) = ^) + C ie ^ (1 - tanh 2 (i)) V ™ , (102) 
where Ci are real constants, z = 1, 2, 3, which satisfy the following relation: 

C 1 + C 2 + C 3 = 0. (103) 



We conclude that exact solutions obtained in [60] are stable in the Bianchi I metric at 
m p < 1/2 and unstable at m 2 > 1/2. The case of m 2 = 1/2 needs a more detailed analysis. 
The first corrections are bounded. 
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6.3 Model with a lump solution 

In the previous subsection kink solutions were considered. In this subsection we consider the 
stability of a lump solution in the model [30] which is motivated by a description of D-brane 
decay within the string field theory framework. We take the one-field cosmological model 
with the potential 

V«) = 2 { l-M>- W- 1 ™ + 3 *r (104) 



and C — — 1. The Friedmann equations (194|) have the following exact solution [30J: 

0o = sech 2 (t), (105a) 

H ^ = 2(3 + 2cosh(t) 2 )tanh 3 (^ 
15m 2 cosh(t) 2 

At t — > oo solution f)105p tends to a fixed point: 

h ' = tL; *' = a (106) 

It is easy to see that 

1^(0) = 0, l&(0) = 4 (l - J^) . (107) 

Using (1551) . we obtain that solution fl 1 5 [) is attractive in the Bianchi I metric at m 2 < 2/5. 
In [30] the authors consider a model without the CDM, at the same time, the results of 
Section 2 show that solution ( 11051) is stable with respect to small fluctuations of the initial 
value of the CDM energy density as well. 

Let us perturb the Friedmann equations in the standard way, 

H = H (t) + eh(t), <f> = </> (t)+e<p(t). (108) 

To first order in e we have the following system of equations: 
2 

h H - sech 2 (t) tanh 2 (t)0 = 0, 

m 2 

\ ( -(4 + cosh(2t))sech 2 (t)tanh 3 (t)/i + (6sech 4 (t) - 4 sech 2 (t))<^ ) - (109) 
m p \5 / 

- 4 ( 2 + 3sech (*)) ( tanh 4( t ) _ 2tanh 6 (t))v9 + 2sech 2 (t) tanh(t)0 = 0. 
25m 2 

System (11091) has the following solutions: 

1 / -4+30m2 / 2cosh 2 (t)-3 \ 

<p = — ^— 5C 2 m 2 cosh(t) 5m 'p j e V 10m p c ° sh4 wy - 2C 1 cosh 2 (t) - 

2 sinh(t) cosh (t) \ 

- 2C 2 I \— (-15m 2 cosh 4 (t) + 10m 2 cosh 6 (t) + 8 cosh 2 (t) - 6 - 4 cosh 6 (t) + 2 cosh 4 (t)) x 

J sinh (i) 
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x cosh(t) ( bm i > e V^^Ht)) ( cos h 2 (t) - m 2 p )dt + 2d 



h _ 16(cosh(2t) - 1) / /"cosh(t) 5m 'p ; e V 10m p cosh4 « 

~ cosh(6t) + 6cosh(4t) + 15 cosh(2t) + 10 V 1+ 2 J sinh 3 (t) X 

x [-15ml cosh 4 (t) + lOmJ cosh 6 (t) + 8 cosh 2 (t) - 6 - 4 cosh 6 (t) + 2 cosh 4 (t)]dt) . 



Using (1105bj) we get 

4 sinh(t) 2 (3+cosh(t) 2 ) 

hi = h+ Ci coshtTf^e 10m p cosh (*' 4 , (110) 

where (7, are arbitrary real constants which satisfy f)103p . 

It is easy to verify that h(t) and hi(t) are bounded functions for any values of the pa- 
rameters. Taking into account that 

/W W -8\ 
*->°° F yi0m 2 p cosh\t) J V ; 

we obtain that <p is bounded at m 2 , ^2/5 and unbounded at m 2 > 2/5. The stability in the 
case of m 2 = 2/5 cannot be analysed without using high order corrections. 

In the examples considered in this and the previous subsections, the potentials depend 
on m 2 , which results in the fact that extrema of potentials are either minima, maxima, or 
inflection points. In the next two examples we consider the opposite case, when potentials 
do not depend on m 2 . 

6.4 Model with a massless phantom field 

At present nonlocal cosmological models are being studied very actively [32], [30] - In 
this and the next subsections we consider the stability of solutions of local models which 
correspond to the nonlocal model with a quadratic potential. These solutions have been pre- 
sented in [36], where the method of localizing the nonlocal model with a quadratic potential 
has been proposed. 

Let us consider the one-field model with zero potential, V(<p) = 0. From the Friedmann 
equations fl94|) we obtain 

3# 2 = ^0 2 + 4 ( 112 ) 
H = (113) 
At A > and C < there are the following real solutions: 



<f> (t) = ± \l - ^arctan ( sinh ( J ^(t -t )]]+ d, 
/ A 

H (t) = 





(114) 
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where to an d C\ are arbitrary real constants. 

Let us consider the stability of the solution (H ,(p ). Substituting H and </> into flBTj) - 
(1621) . to first order in e we obtain 



<f(t) = ± 1 , — == r C 3 + C 2 , 



h(t) 



where C 2 and C3 are arbitrary real constants. It is obvious, that functions h(t) and (pit) are 
bounded. In the Bianchi I metric we have 



(115) 

2C 3 



cosh (24P%- to) 1 + 1 



hi{t) = hit) + Ci 



I tanh 2 [ V^^(t_to) I. (116) 



T v 

where real constants Ci, i — 1, 2, 3, satisfy the following relation: 

Ci + C 2 + C 3 = 0. (117) 

Thus, we have obtained that the kink-type solutions (11141) in the Bianchi I metric have 
the bounded first corrections. 



6.5 Model with a quadratic potential and the cosmological con- 
stant 

Let us consider the model of a scalar field with a quadratic potential and the cosmological 
constant. In this case the Friedmann equations are 

H , = ^,^ 2 + B^ 2 + A y (n8) 

H= -47rG N C<j) 2 , (119) 

where C and B are arbitrary nonzero real numbers. 

System (I118l) - (lll9p has the following particular solutions: 

H (t) = ht, Mt) = ht, (120) 

where 



B %2 B 

•2 



*>=-#•< **=12^' (12i; 
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From (I12ip it follows that the function is real if and only if B > 0. The above-mentioned 
solutions exist only if 

A =~^kc- (122) 

To analyse the stability of these exact solutions, we substitute 

H(t) =k x t + £h(t) (123) 

and 

<j){t) = k 2 t + £ip{t). (124) 

in (HHD and ffTT9|) . 

To first order in e we obtain the following system of equations: 

m = - f (*«> + Sg^5*<«)) • (125) 

h(t) = -8irG N Ck 2 ip(t). (126) 

Solutions of (H2SD— (H2SD are 

h{t) = D ie £ t2 + D 2 , (127) 

^) = " a + ^e^* 2 ) , (128) 

where D\ and -D2 are arbitrary constants, 



* =± Vi2sb* (129) 

Therefore, the functions h(t) and <p(t) are bounded at C/B < 0. Real solutions exist 
only if I? > 0, and hence, C < 0. We come to the conclusion that solution (11201) can be 
stable (the first corrections are bounded) only if C < 0, in other words, <p(t) is a phantom 
scalar field. In this case H (t) > at t > 0, hence hi are bounded as well. Indeed, 

= + C ie ' 3 o H ° {T)dT = (A + a)e^* 2 + £> 2 , (130) 
are bounded. Constants Cj satisfy the relation (I103p . 

7 Conclusion 

We have analysed the stability of isotropic solutions for the models with NEC violation in 
the Bianchi I metric. 

In our paper for the one-field model with the CDM we used the Lyapunov theorem and 
found sufficient conditions for stability of kink-type and lump-type solutions both in the 
FRW metric and in the Bianchi I metric. The obtained results allow us to prove that the 
exact solutions, found in string inspired phantom models [6TJII3T)] . are stable. A generalization 
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of this result to two-field models, for example, quintom models, requires further studies and 
will be considered in future investigations. 

We found the explicit form of the connection between hi(t), h 2 (t), and h 3 (t), which define 
metric perturbations in the Bianchi I metric, and ho, which defines perturbations in the FRW 
metric. We have proved that fluctuations for the fields and the CDM energy density in both 
metric are the same. In particular, for Hq ^ the boundedness of ho is a sufficient and 
necessary condition for the boundedness of hi(t), h 2 (t), and h 3 (t). This result is valid for 
both iV-field and /c-essence models. 

Note that linear-in-time solutions for the simple models with a quadratic potential, which 
have been considered in Subsection 6.5, are stable with respect to the first corrections only 
in the phantom case (which corresponds to H > 0). It means that NEC violation does 
not lead to the instability in this sense. It gives us an intuitive reason to expect that more 
complicated NEC violated models can also have stable isotropic solutions. This expectation 
has been confirmed in the models considered in Section 6. 

To conclude, our results are for the NEC violated models. If energy conditions are 
satisfied, then results, similar to those obtained in Section 4, are consequences of the Wald 
theorem [55] and its generalizations [56J, [57J EH [59] . 

Our study of the stability of isotropic solutions for the models with NEC violation in the 
Bianchi I metric shows that the NEC is not a necessary condition for classical stability of 
isotropic solutions. Because of strong limits on anisotropic models from observations [631164]. 
cosmological models developing large anisotropy should be discarded. In this paper we have 
shown that the models [301 ESS ED] have stable isotropic solutions and that large anisotropy 
does not appear in these models. 
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